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Abstract 

We study various aspects of the topological quantum computation scheme based on the non- 
Abelian anyons corresponding to fractional quantum hall effect states at filling fraction 5/2 using 
the Temperley-Lieb recoupling theory. Unitary braiding matrices are obtained by a normalization 
of the degenerate ground states of a system of anyons, which is equivalent to a modification of the 
definition of the 3- vertices in the Temperley-Lieb recoupling theory as proposed by Kauffman and 
Lomonaco. With the braid matrices available, we discuss the problems of encoding of qubit states 
and construction of quantum gates from the elementary braiding operation matrices for the Ising 
anyons model. In the encoding scheme where 2 qubits are represented by 8 Ising anyons, we give 
an alternative proof of the no-entanglement theorem given by Bravyi and compare it to the case of 
Fibonacci anyons model. In the encoding scheme where 2 qubits are represented by 6 Ising anyons, 
we construct a set of quantum gates which is equivalent to the construction of Georgiev. 
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I. INTRODUCTION 



Quantum computers are expected to be able to perform calculations which are impos- 
sible for classical computers, due to quantum entanglement and quantum parallelism 
Unfortunately, quantum computers seem to be extremely difficult to build because of the 
unavoidable noise and decoherence caused by the coupling of the qubits and the ambient en- 
vironment. It is a daunting task to construct a quantum computer that has a large number 
of qubits and has a low error rate. However, there is a promising approach, called topologi- 
cal quantum computation (TQC) |2-|l3|. proposing to encode the qubit information into a 
topological quantum field. Kitaev [2] proposed that a system of anyons can be considered 
to be a quantum computer. Unitary matrices are related to moving the anyons around each 
other. Measurements are performed by joining anyons in pairs and observing the result of 
fusion. Interference experiments 14H19| are also proposed to initialize and read out quantum 



states. The computation is fault-tolerant by the topological nature of the system. 

Different from fermions and bosons, which are the totally antisymmetric and the sym- 
metric representations of the permutation group S n , anyons carry fractional charges and 



have fractional statistics 



20 



2f| which result in nontrivial phases (for Abelian anyons) or 



matrices (for non- Abelian anyons) for permutations. In fact, the underlying symmetry of 
the system of anyons is the braid group B n . Abelian anyons correspond to one-dimensional 
representations of B n and the quantum gates one can construct from them are very limited 
221 ] . Non- Abelian anyons, on the contrary, are much more useful to TQC, since the braiding 
of non-Abelian anyons induces non-commuting (non-Abelian) representations of B n , from 
which one can construct various quantum gates. 

Physically, anyons are collective excitations in some condensed matter systems, such as 
the fractional quantum hall effect (FQHE) states of two dimensional electron liquids. For 
example, the effective theories of FQHE states with filling levels v = 1/3, v = 5/2 and 
v = 12/5 correspond to Abelian anyons, non-Abelian Moore- Read 23] and Read-Rezayi 



241 ] anyons respectively. Mathematical 



Chern-Simons effective field theories 



y, properties of anyons can be described by SU(2)k 



251 ] and parafermion conformal field theories (CFT) 



261 ] . the k = 2 and k = 3 cases corresponding to the v = 5/2 and the v = 12/5 FQHE states 



respectively. 

The SU(2)k Chern-Simons theory is a topological quantum field theory [25J] which has 
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a deep relationship to knot invariants, Jones Polynomial [27] especially. Kauffman and 



Lomonaco 



28 



29| studied unitary representations of braid groups in terms of g-deformed 



spin networks, or Temperley-Lieb recoupling theory 30]. 

It is one of our purpose to apply the method of Kauffman and Lomonaco to calculate 
explicitly the elementary braiding operation (EBO) matrices which govern the exchanges 
of Ising anyons, the first non-Abelian anyons model proposed by Moore and Read 23] by 
constructing a wave function (the Pfaffian state) for the v — 5/2 FQHE state corresponding 
to the 577(2)9. Chern-Simons theory. Direct experimental observation of fractional electron 
charge e/4 31, 32]] at the v = 5/2 FQHE state gives some evidences in support of the 
non-Abelian nature of this state. There are many works concerning the braiding properties 



of the Ising anyons using CFT method [33j, |34| or quantum group method 3^. Nayak and 



Wilczek [33] suggested that the Pfaffian wave functions of n (n even) Ising anyons form a 
2 n / 2_1 dimensional spinor irreducible representation of the rotation group SO(n), to which 
a rigorous treatment is given by Georgiev recently 34]. Quantum group approach 35] also 
gives equivalent results. As we will see, the EBO matrices for the Ising anyons can be 
elegantly derived by using the Temperley-Lieb recoupling theory. 

One of the attractive properties of the Ising anyons TQC model is that the excitation 



gap at the corresponding filling fraction v = 5/2 is the highest one among a 



1 non-Abelian 



FQHE states, resulting in a very low (or even lower) error rate of 10 30 36H38]. Although 
this Ising anyons model is not universal fl-EJ for TQC, i.e., the braid group representations 
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45|. 



are not dense in unitary groups, it receives extensive attention in the past few years 
In fact, it is proved by Bravyi |39| that no entangled states in the computational space can 
be obtained purely topologically and the Ising anyons TQC model is classically simulatable. 
We show that the same conclusion can be obtained from the Temperley-Lieb recoupling 
theory approach. 

We should stress that this no-entanglement theorem does not mean that there is no 
entanglement between Ising anyons at all. This rule only applies to the qubit encoding 
scheme (which is consistent with the quantum circuit model) where each qubit is encoded 
in 4 Ising anyons. Entangled quantum gates, such as controlled-Z and controlled-NOT 
(CNOT), can be realized purely by braiding 6 Ising anyons with definite topological charge 
(quantum spin). Using the EBO matrices obtained in this paper, we construct a set of 
useful 1-qubit and 2-qubit quantum gates, which is not the same as, but equivalent to the 
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construction by Georgiev 4l|, |42| . 



The outline of the remainder of this paper is as follows. Section 2 reviews the general 
models of non-Abelian anyons for TQC and the formalism of the Temper ley- Lieb recoupling 
theory needed later. Unitary representations of the braid group in the Hilbert space of the 
degenerate ground states of non-Abelian anyons are obtained by a physical argument which 
requires that the fusion paths of the anyons form an orthonomal basis of the Hilbert space. 
In section 3, we derive the EBO matrices of the Ising anyons model. In section 4, we study 
some aspects of Ising anyons TQC using the results of section 3. Conclusions and discussions 
are presented in Section 5. 

II. TEMPERLEY-LIEB RECOUPLING THEORY AND UNITARY REPRESEN- 
TATIONS OF BRAID GROUPS 

In this section, we first review the definition of the quantum states of a system of anyons 



and the Temperley-Lieb recoupling theory [30| and then discuss the method to produce 
unitary representations of braid groups. 

A. Models of non-Abelian anyons for TQC 



A model of non-Abelian anyons consists of the following three elements [llj: a list of 
particle types, the fusion rules, and the braiding rules. In the formalism of SU(2)k Chern- 
Simons theory, anyons are quasi-particles having half-integer g-spins (spins for short) s = 
0, 1/2, 1, • • • , k/2 as their quantum numbers. The fusion rules of these particles are truncated 
versions of the rules of addition of ordinary angular momenta, 

si <g> s 2 — |si — S2I ffi \sx — s 2 \ + 1 © • • • © min(si + s 2 , k — s x — s 2 ). (1) 

When a number of non-Abelian anyons with definite spins fuse consecutively into a single 
anyon with some spin, the sequences of the intermediate spins of the fusion paths represent 
different quantum states of the Hilbert space. 

Anyons commonly appear as collective excitations in 2 dimensional systems. When they 
move, their world lines propagate in a 3 dimensional space-time. Thus the exchange of a 
pair of anyons corresponds to the braiding of their world-lines. (We will call the braiding 
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of the world-lines of anyons shortly as the braiding of anyons, but it is important to keep 
in mind what it means actually.) In TQC, we perform quantum computations by braiding 
anyons to realize certain quantum gates. Any braiding can be expressed as a sequence of 
EBOs whose representations in the above Hilbert space are the EBO matrices we want to 
find. The essential task of deriving the EBO matrices is the determination of the so called 
R-matrix and F- matrix introduced first in the context of CFT [49] . The former is the unitary 
matrix inducing the exchange of neighboring anyons with definite total spin, and the latter 
accounts for the associativity of fusions of anyons. In the next subsection, we will give their 
diagrammatic definitions in terms of the Temperley-Lieb recoupling theory. 

B. Temperley-Lieb recoupling theory 



Temperley-Lieb recoupling theory [30J is based on the Kauffman bracket polynomial 
model for the Jones polynomial at roots of unity and the tangle-theoretic Temperley-Lieb 
algebra. 

1. Braid group and Temperley-Lieb algebra 

The Artin braid group B n can be presented as a set of generators 01, 02, o n -\ that 
obey the following relations, 



The Temperley-Lieb algebra TL n can be presented similarly as a set of generators 
Ui, U 2 , U n _i, whose representations are related to the representations p(o"j) of B n by 



where the Kauffman variable A is taken to be A = ie ln ^ 2r for Jones polynomial at 4r-th 
roots of unity such that the quantum dimension of the spin 1/2 anyon is d = —A 2 — A~ 2 = 
2 cos(7r/r). 




OiOj = OjOi for \i — j\ > 2; 

= for i — 1,2, 



n-2. 



(2) 



p{a l ) = A + A- 1 U l 



(3) 



5 



c 



FIG. 1: The Jones-Wenzl projector and the 3-vertex. The integers i, j, and k are determined by 
the relations a = j + k,b = k + i, and c = i + j. 

2. R-matrix and F-matrix 

The basic object of the Temperley-Lieb recoupling theory is the Jones-Wenzl projector 



301 ] . The left graph in Fig. [T] shows the Jones-Wenzl projector constructed on the basis 
of the Kauffman bracket polynomial expansion. The n-strand projector corresponds to the 
world line of an anyon with spin n/2. The middle graph in Fig. [1] shows the 3-vertex 
constructed from the projectors which corresponds to the interaction (fusing or splitting) of 
3 anyons with spins a/2, 6/2, and c/2. The right graph in Fig. [1] is a simplified notation 



for the 3-vertex. Note that the g-admissible conditions 



30] for the 3-vertex, 



a + b + c = even; 

a + b-c>0, b + c-a>0, c + a-b>0; (4) 
a + b + c<2r~A 

say exactly the same thing as the fusion rules Eq. ([!]) do due to the relation r = k + 2 between 



SU(2)k Chern-Simon theory at level k and Jones polynomial at 4r-th roots of unity [25] and 
the fact that the projector with label n represents the world line of an anyon with spin n/2. 
Various spin networks can be constructed from Jones-Wenzl projectors and 3-vertices. See 
Fig. El 

As the world lines of anyons and their fusions being identified with projectors and 3- 
vertices, we now consider the braiding properties of anyons in the context of the Temperley- 
Lieb recoupling theory. Fig. [3] shows the braiding of two anyons with spins a/2 and b/2 



A 



n 




Tfa b i\ 

1 \cdj) — 



n 





6(a, b, c) 





8(a,6,c) r 



C 



C 



FIG. 2: The definitions of the delta net A n , the theta net Q(a,b,c), and the tetrahedron net 



r(*3}). Formulae 



30f | for evaluating these spin-nets are presented in the appendix. 
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FIG. 3: The action of the R-matrix, whose matrix elements are given by Eq. ([5]). Our convention 
is that seen from above, the braiding of a and b is counterclockwise. 

fusing into a spin c/2 anyon. Since this operation dose not change the total spin of the two 
fusing anyons, the corresponding matrix, the R-matrix is diagonal in the underlying Hilbert 
space. The matrix element R^ b is given by the following formula jsoj], 



r: 



ab 



/_-^\(a+b~c)/2^-(a(a+2)+b(b+2)-c(c+2))/2 



(5) 



Not all braids are of this case, in which the two braiding anyons fuse into a single anyon. 
To see this, it is sufficient to consider the case of 4 anyons with total spin 0, which is 
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FIG. 4: The configuration of 4 anyons a, b, c, and d fusing into (the vacuum) is equivalent to the 
configuration of 3 anyons a, b, and c fusing into one anyon d. The first equality follows from the 
fact that one can get a only by fusing an anyon with label j = d with the anyon with label d. 
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FIG. 5: The use of the F- matrix. 



equivalent to 3 anyons fusing into the 4-th anyon. This is shown in Fig. HI 

The braiding of b and c in Fig. [5] can not be accomplished via a single R-matrix and 
is realized only by a combination of the R-matrix and the F-matrix. Fig. [6j shows the 
definition 30(| of the matrix element F(^)ij of the F-matrix as well as the formula to calculate 
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= FCZ) 



\cd)ij 





b KcdJij 




z?(ab\ S{a,b,k)e(c,d,k) 
r KcdJik 



A; 



FIG. 6: The upper diagram shows the definition of the matrix element F(cd)ij of the F-matrix. A 
summation on j is implied. The lower diagram shows the derivation of the expression Eq. ([6]) of 
the F-matrix in terms of the delta net, the theta net, and the tetrahedron net. Expressions in Fig. 
[2] are used in the derivation. 



it in terms of the delta net A n , the theta net 9 (a, b, c), and the tetrahedron net T{^\ 1 -) 



©(a, 6, j)Q(c, d,j) 



(6) 



C. Unitary representations of Artin braid groups 

However, the F-matrix denned above in the Temperley-Lieb recoupling theory is not 
unitary, resulting in a non-unitary representation of the braid group. Unitary (in fact real 
and orthogonal) F-matrix, and hence unitary representation of B n can be obtained by a 
redefinition for the basis states in the Hilbert space of the anyons. The guideline of the 
following argument is the requirement that the fusion paths should represent an orthonormal 
basis of the Hilbert space of the degenerate ground states of a system of anyons. 

We need only to consider the orthonormal problem of the states of four anyons fusing 
into the vacuum. The definition of the state \i) and the calculation of the inner product 
of the two states \i) and \j) are shown in Fig. [7J We see that the orthogonal property is 
already satisfied oc 6^), but the state vectors are not normalized ^ 1). It follows 
that an orthonormal basis {\i}} = <%) of the Hilbert space is obtained by normalizing 

9 




FIG. 7: The upper diagram shows the graphic definition of the state vector \i) for 4 anyons with 
total spin 0. The graphic representation of the dual vector (i\ is defined to be the upside-down 
version of that of the vector The lower diagram shows the calculation of inner product of 
two states \i) and \ j). Expressions in Fig. [2] are needed in the calculation. 
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b d 



i) = 



S(a,c,i)S(b,d,i) 




FIG. 8: The normalized basis {\i)} is obtained from the original one {\i}} by multiplying an 
appropriate factor to each state 

each of the states, as depicted in Fig. |HJ In the orthonormal basis {\i)}, the new F-matrix 
can be derived to be 



rpia bi\ 
zl \cdjj 



(7) 



'©(a, b, j)©(c, d,j)Q(a, c, i)Q(b, d, i) 

This new F-matrix is real and orthogonal (hence unitary), as we will see in explicit calcula- 
tions latter. We note that Kauffman and Lomonaco [28] obtained the unitary F-matrix by 
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FIG. 9: The fusion diagram for n + 3 (n odd) Ising anyons with total spin 0. The first (the left 
most one) and the second Ising anyons fuse into an anyon with spin x\ = or 1, and then x\ fuse 
with the third Ising anyon into an Ising anyon with spin 1/2, • • -, and finally, x n fuse with the 
(n + 2)-th Ising anyon into the (n + 3)-th Ising anyon. 

multiplying each 3-vertex with the following factor, 

/(a, b, c) = (A a A b A c ) 1 /4 /(0(G) 6j c)) i/2_ (g) 

One can check that this modification to 3-vertices results in the same F-matrix given by 
Eq. ©. 



III. EBO MATRICES FOR ISING ANYONS 



In this section, we apply the results of section 2 to obtain explicitly the representations of 
the generators of the braid group governing the exchanges of Ising anyons. Our presentation 



28 



29]. 



follows closely to Kauffman and Lomonaco 

For the Ising anyons model, the Kauffman variable is A = ie ln ^ 8 , and the quantum 
dimension of the spin 1/2 Ising anyons is d = -A 2 -A- 2 = y/2. The allowed spins of anyons 
in this model are 0, 1/2 and 1 and the fusion rules for these anyons can be deduced from 
Eq. CQ) to be, 

0®j=j for j =0,1/2,1; 

1/2 ® 1/2 = 1; (9) 
l®j = l-j for j = 0,1/2,1. 
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The dimension of the Hilbert space of n spin 1/2 Ising anyons with total spin is 2™/ 2 ~ 1 , 
approaching d n in the limit n — > oo. According to the above fusion rules, the fusion diagram 
of n + 3 (n must be odd) Ising anyons with total spin takes the form as shown in Fig. [9] 



50] 



We now calculate the unitary representation of B n+3 for the braiding of the n + 3 Ising 
anyons. We denote the elementary braiding operation (EBO) of the first and the second 
Ising anyons as a\, the EBO of the second and the third Ising anyons as a 2 , The 
corresponding EBO matrices are denoted by p[pi) where 1 < i < n + 2. 

The first EBO matrix p(<Ji) is easy to calculate. It depends only on the label x\. In the 
basis {x\ = |0), |1)}, p{(J\) is simply given by the following R-matrix, 

/< o j J-** A (10) 
[on?) { o a) 

which corresponds to the following Temperley-Leib generator, 

U = AR-A 2 =\ = • (11) 

v°°y I °) 

To calculate p(02), we need to calculate the following F-matrix in the same basis as above 
(using Eq. ([7]) and the formulae in the appendix), 

F = ( F{lllU F{ll) ° 2 )=±( 1 1 ) (12) 
\F(\\U F(\\) 22 ) V2[l-l)' 

corresponding to the following Temperley-Leib generator, 

'-(£#-*(::)■ 

where S = p{a 2 ) = FRF- 1 . 

Now consider the case of p(cr 3 ). When either x\ or x% or both of them equal to 0, the 
situation is similar to the case of p(<Ti) where only R-matrix elements are needed to be 
calculated. The case in which both xi and x% are 1 deserves special consideration. The 
EBO 03 does not change the value of x 2 when both x± and x% are 1, and the matrix element 
of p((Jz) in this case is found to be the same as in the case where both x\ and x% are by 
doing some graphic calculations. 
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The other EBO matrices can be calculated in the same way. By choosing the basis of the 
n + 3 Ising anyons with total spin (n odd) as {\x\, x 2 , ■ ■■x n )} where Xi equals to 0, 1/2, or 
1 such that the fusion rules are met at each fusion vertex along the whole fusion path, as 
shown in Fig. [HI we can find a representation of the Temperley-Lieb algebra TL n+3 . For the 
first and the last generators of TL n+3 , we have, 

U 1 \0,x 2 ,...,x n ) = y/2\Q,x 2 , ■~,x n ); 
Ui\l,x 2 , -,x n ) = 0; 

(14) 

U n+2 \x Xl ...,x n _i,0) = \/2|xi, ...,x n _i,0); 
1) = 0. 

For U 2 and U n+ i, we have, 

U 2 \0,l/2,x 3 ,...,x n ) = 

75IO, 1/2,23, -,Xn) + 1/2,23, ...,2 n ); 

C/ 2 |l, 1/2,23, -;X n ) = 

75IO, l/2,2 3 ,...,2 n ) + ^|1, 1/2,23, ...,2 n ); ^ 
C/ n+ i|2i, ...,2 n _ 2 , 1/2,0) = 

^|2i, ...,2 n _ 2 , 1/2,0) + ^|2i, ...,2 n _ 2 , 1/2, 1); 
f/ n+ l|2i, ...,2 n _ 2 , 1/2, 1) = 

^\xx, ...,2 n _ 2 , 1/2,0) + ^\X1, ...,2 n _ 2 , 1/2, 1). 

For the middle ones, we have, 



Ui\xi, .. 


• 5 Xi- 


.3,0, l/2,0,2 i+ i, .... 


Xn) 


V2\x u . 


. . , 2j 


-3,0, l/2,0,2 i+ i, .. 




Ui\x u .. 


., 2j_ 


-3) 0, 1/2, 1, 2 i+ i, .... 


x n ) = 0; 


Ui\xx, .. 


• ; Xi- 


.3, 1/2,0, 1/2, 2 J+1 , 


• • • , ^n) 




■ ■ , 2j. 


_ 3 , 1/2,0, 1/2, 2 i+1 , 


• • • , Xfij ~\~ 




• • , 2j- 


_ 3 , 1/2,1, 1/2, x i+ i, 


2 n ) , 


t/"i|2i, .. 


• 5 2j_ 


.3, 1/2,1, 1/2, x i+1 , 


• • • , 2 n ) 




■ • j 2j. 


_ 3 , 1/2,0, 1/2, x i+1 , 


2 n )-|- 




• • ) 2j- 


_ 3 , 1/2,1, 1/2, 2 i+1 , 


2 n ), 


f/i|2i, .. 


., 2j_ 


-3, 1, 1/2, 0, 2 i+ i, .... 


2n) = 0; 


Ui\x ± , .. 


• 5 Xi- 


-3) 1) 1/2, 1, 2j + i, .... 


Xn) 




.., 2j 


-3, 1, 1/2, 1, 2 i+ i, .. 


■ , x n ) ■ 



(16) 
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The EBO matrices can be obtained immediately from the above representations of the 
Temperley-Lieb algebra by using Eq. ([3]). One can check that these representations of the 
braid group indeed satisfy the Artin relations in Eq. (j2J). 



IV. CONSTRUCT QUANTUM GATES FROM THE EBO MATRICES FOR THE 
ISING ANYONS MODEL 

After obtaining the EBO matrices of the Ising anyons model, we study in this section 
some aspects of Ising anyons TQC. 

The qubit encoding scheme which is consistent to the quantum circuit model is to use 
each group of 4 Ising anyons with total spin for each qubit such that an n-qubit system 
uses 4n Ising anyons. This is the encoding scheme used by Bravyi [39| who proved a no- 
entanglement theorem which states that entangled 2-qubit states can never be prepared by 
pure topological braiding operations. The proof by Bravyi uses the stabilizer constrains and 
the no-leakage error conditions. In the following, we give a graphical demonstration of this 
result from the Temperley-Lieb recoupling approach. 

In this qubit encoding scheme, 2-qubit states are encoded in 8 Ising anyons, 4 Ising anyons 
for each qubit. See Fig. [HH The two groups of EBOs {a±, a 2 } and {oe, 07} apply completely 
within the first and the second qubits respectively and can not generate entanglement be- 
tween the two qubits. Since ex 3 depends on x\ and x 3 , a 5 depends on x 3 and x 5 , and cr 4 can 
change x 3 from to a superposition of and 1, it is possible to create an entangled state by 
a sequence of these EBOs, such as a% l a± a 5 a^a 3 . However, it is impossible to avoid leakage 
errors by braiding this way. To see this, it is convenient to change the fusion paths in Fig. [10] 
to another basis as shown in Fig. [TT], where the total spin of the 2-qubit system can be either 
or 1. The EBO sequence 03 crj a^a^az in Fig. [TD] is equivalent to the single EBO cr 3 in 
Fig. EJ Since the braid matrices of 03 in Fig. [11] for the two sectors (total spin and 1) are 
not equivalent, and entanglement can not be created without using cr 3 in Fig. [11] leakage 
error from the computational space (labeled by x\ and x§ in both Fig. [10] and Fig. [TT]) to 
the uncomputational space (labeled by 23 in Fig. [10] and x to tai in Fig. [TT]) is unavoidable. 
It is instructive to compare the above situation with the case of the Fibonacci anyons 



ll| model which is universal for TQC 4j-[6j]. See Fig. [10] too. For Fibonacci anyons, each 



intermediate spin Xj can be either or 1 (as long as no two 0s appear consecutively) and 
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(a)abcdefgh 

w w w w 



CTi CT 2 CT 3 (7 4 (J5 CT 6 (T7 

(b) a 6 c d e f 9 h 



X\ X 2 Xs X4 x 5 

FIG. 10: Part (a) shows the initial state of the two qubits. The 4 anyons a, b, c, and d form the 
first qubit and the other 4 anyons e, /, g, and h form the second qubit. Part (b) shows the full 
quantum labels needed when we braid the world lines of the anyons. The computational space is 
spanned by the states { \x\,x§) = 1 0, 0) , |0, 1), |1,0), |1,1) }• In this notation, the initial state in 
part (a) is 1 0, 0) . For Ising anyons, 22 and X4 can only be 1/2; for Fibonacci anyons, X{ can be 
either or 1 such that x% + > 0. 



a b c j g h a b c j g h 




FIG. 11: The formation of the spin (x to tai = 0) and spin 1 (a^otai = 1) sectors for 8 Ising anyons. 
Each sector is 4 dimensional and a basis { \x\,x£) = |0, 0), |0, 1), |1,0), |1,1) } can be chosen for 
both sectors. The rightmost diagram in this figure is related to part (b) of Fig. [10] by a change of 
basis using a F-matrix. 
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FIG. 12: The single-qubit states with 4 Ising anyons and the braiding diagram for the Hadamard 
gate. x\ = and x\=l correspond to the states |0) and |1), respectively. 
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entangled states can be generated by braiding the 8 Fibonacci anyons. Bonesteel et al 
481 ] constructed some entangled 2-qubit gates such as the controlled-iX gate by weaving two 
Fibonacci anyons (c and d) from the control qubit into the target qubit which approximates 
the identity matrix, followed by a braiding within the target qubit which approximates the 
iX gate, and then weaving them back to their original positions. The controlling operation 
is realized by virtue of the fact that the whole braiding does nothing when the total spin of 
c and d is and acts as the iX gate on the target qubit when the total spin of c and d is 1. 
The nonuniversality of the Ising anyons model prevents us from realizing entangled gate in 
this way. 

We note that the no-entanglement theorem only applies to the above qubit encoding 
scheme where 2 qubits are represented by 8 Ising anyons. Entangled quantum gates can be 



constructed in a different qubit encoding scheme, as studied by Georgiev [41|, |42|. In this 
scheme, 1-qubit and 2-qubit states are encoded in 4 and 6 Ising anyons with total spin 
respectively. 

Consider 1-qubit gates first. Taking the basis of the Hilbert space to be (see Fig. ^ 
{\xi) = |0), |1)}, the two dimensional EBO matrices of for four Ising anyons are found 
to be 



p (2 Vi) = p (2 V 3 ) 

—i-w/Si ( 



'-10 



p (2 V 2 ) 



V 

1 z 
i 1 



% 



(17) 



(18) 



V2 

One can construct the Hadamard gate H, the phase gate S, and the three Pauli gates 
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X, Y, and Z using the two dimensional EBO matrices given above (~ means equal up to 
an unimportant global phase), 



H 




p {2 \(T 1 a 2 a 1 ) 



(19) 



s 



X 



1 

i 

1 

1 



Y 




P (2 Vr x ); 



p [2 \o 2 a 2 ) 



p (2 Wi)- 



(20) 



(21) 



(22) 



(23) 



Fig. H2] shows the encoding of the 1-qubit states as well as the Hadamard gate constructed 
by three braids. 

However, it fails to construct the 7r/8 gate, 

^ I U 1 If U I 

(24) 



T 








,,471-/8 



reflecting the fact that the Ising anyons model is not universal for quantum computation. To 



39. 



40]. 



remedy this, we have to supplement braiding with some non-topological operations 

Now consider the 2-qubit case. Taking the basis of the Hilbert space to be (see Fig. [H]) 
{ \x 1 ,x 2 ,x 3 ) = |0, 1/2,0), |0, 1/2,1), |1, 1/2,0), |1, 1/2,1) }, the four dimensional EBO 
matrices of Bq for 6 Ising anyons with total spin read 



p( 4 V0 = e-/ 8 diag(-l,-l,M); 



( I i 0^ 
010i 
i 1 
Oi 01 



P {i) M 



V2 



,(4), 



a 3 ) = e^ /8 diag (-l,i,i,-l) 



(25) 



(26) 



(27) 
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P (4 V 4 ) 



V2 



( I i 0^ 
l 100 

001 ! 

% 1 
-l,z, -l,z) . 



(28) 



/ 



e <ir/8 diag(-l,i,-l,i). (29) 
From the four dimensional EBO matrices given above, one can construct useful 2-qubit 



quantum gates, such as CNOT (up to an unimportant global phase), 

^ 1 ^ 



CNOT 



10 
1 
10 



p (4) (cr 3 V 4 V 5 1 cr 3 o-4a 3 cri) 



(30) 



and controlled- Z 



^ 1 ^ 



controlled- Z 



\ 



(31) 



/ 



10 
10 
-1 

Fig. [13] shows the encoding of the 2-qubit states as well as the braiding diagram of 
CNOT. Note that the braid sequences for CNOT and controlled- Z are not unique, which 
is a consequence of the Artin relations Eq. ([2]) for the generators of the braid group. Our 



41 



42]. 



constructions are different from, but equivalent to the ones given by Georgiev 

Note that the 4 dimensional Hilbert space of 6 Ising anyons with total spin is a subspace 
of the 8 dimensional space of 8 Ising anyons. The other subspace, which is also 4 dimensional, 
corresponds to 6 anyons with total spin 1. See Fig. [TT] Using the method in section 3, we can 
also find the EBO matrices for the spin 1 sector. It turns out that the EBO matrices p^'(ai), 
P^( a 2), p^(cr^), and p^(cr 5 ) in the spin 1 sector take the same form as in the spin sector, 
but p^((Js) has a different form, p^(o" 3 ) = e l?r//8 diag (i, — 1, — 1, i). Therefore, entangled 
2-qubit quantum gates will have different braid sequences in the spin 1 sector. For example, 
one possible braid sequence in the spin 1 sector for CNOT is a^ 1 cr^ 1 aj 1 a^a^ai, which is 
not topologically equivalent to the one given in Eq. ( |30|) . Despite of the non-equivalence of 
the braid sequences for a given quantum gate in these two sectors, the computational power 



of the two sectors are equivalent 



M\ 
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x 3 




FIG. 13: The 2-qubit states with 6 Ising anyons with total spin and the braiding diagram for the 
CNOT gate. x\ and 23 can be or 1 and X2 can only be 1/2. Note that the braid sequence for 
CNOT is not unique, a consequence of the Artin relations for the generators of the braid group. 

V. CONCLUSION AND DISCUSSION 

As demonstrated in previous sections, the Temperley-Lieb recoupling theory provides a 
natural language for describing the braiding properties of non-Abelian anyons. We have 
applied this theory to derive the EBO matrices of the Ising anyons model. We paid a special 
attention to the normalization of the degenerate ground states corresponding to the fusion 
paths of the anyons. This normalization results in the correct unitary F-matrices and is 
equivalent to the redefinition of the 3- vertices proposed by Kauffman and Lomonaco 28]. 

One important feature for the construction of the two-qubit gates is that we can not 
construct them without the use of 03, the EBO acting between the two qubits. This is 
because that the EBOs o\ and 02 act only on the first qubit and 04 and 05 act only on the 
second qubit. Indeed, the first two and the last two EBO matrices can be expressed as a 
tensor product of two matrices, and the middle EBO matrix p^(as) can not, reflecting the 
(topological) entanglement of the 2 qubits. This entanglement is crucial for the construction 
of the 2-qubit entangled gates. However, to get this entanglement, we need to project the B 8 
representation to either the spin or the spin 1 Bq representations. Alternatively, entangled 
quantum gates can be constructed by parity measurement as well as braiding operations 
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39|,|43|. 

The construction of the 2-qubit gates in each sector can be easily achieved by brute 
force search, since the braid lengths of controlled-Z and CNOT are very short (3 and 7 
respectively). However, there is a more heuristic approach, namely, the genetic algorithm 
(GA) approach. A possible braid sequence for the CNOT gate can be found within a minute 
using GA, while it takes a much longer time using the brute force approach. The superiority 
of GA over brute force search is not significant for Ising anyons TQC, but we expect that 
here is a potential application of GA to Fibonacci anyons topological quantum compiling 



461448]- 



Acknowledgements 

We thank Jens Fjelstad and Ben Goertzel for numerous discussions. We also thank the 
referees of EPJB who pointed out a mistake of our original manuscript and helped improve 
this paper a lot. Z. Fan is supported by National Natural Science Foundation of China 
under grant numbers 10535010, 10675090, 10775068, and 10735010. 



Appendix A: Formulae for evaluating the spin-nets 

In this appendix, we present the formulae for the evaluations of the A-net, the G-net, 



and the tetrahedral net 



3Q|. 



The A-net evaluation is 

A n = (-1)> + 1], (Al) 

where [n] is the g-deformed integer defined as [n] = (A 2n — A~ 2n )/(A 2 — A~ 2 ). The 0-net 
evaluation is 

B(a b c) - ( lV W* + .?+fr+l]'ffi]'W! (A2) 

° (a ' 6 ' C) - ( 1} [ i + j ][[ j + k ]\[k + i]\ ' (A2) 

where the g-deformed fractional [n]\ is defined as [n]\ = [n][n — 1]...[2][1], and the integers i, 
j, and k are determined by the relations a = i+ j,b = j + k, and c = k + i. The bracket 
evaluation of the tetrahedral net is 

T (abi\ _ I\ m ,J b "- a ™V v (-l) s [s+l]! /Ao\ 

1 \cdj) [a]\[b]\[c]\[d}\[i]\[j}\ ^max{a m }< s <mm{M Y[ m [s-a m ]\Y[ n [b n -s]\i 
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where a m and b n are given by a\ = (a + d + i)/2, a<i = {b + c + i)/2, a 3 = (a + b + j)/2, 
a 4 = (c+d + j)/2, h = (b + d + i+j)/2, b 2 = (a + c + i+j)/2, and b 3 = (a + b + c + d)/2. 
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